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FOR REACTION-WHEEL ATTITUDE CONTROL* 


** 

Robert H. Cannon, Jr. 


ABSTRACT 


In many space vehicles, attitude control is best accomplished with 
combination systems using reaction wheels for momentum exchange and 
storage, plus jets for periodic momentum expulsion. Design of the 
reaction-wheel control involves evaluating the time history of system 
response to disturbances, many of which are either sinusoidal or impul- 
sive . 


As an aid to such evaluation, this paper developes basic response 
relations --vehicle attitude, control torque, wheel motion, mechanical 
power, and energy consumption — for a vehicle subjected to both types of 
disturbance. Limiting values are calculated, assuming no standby losses. 
(The possibility of exchanging momentiam with minimum energy loss is dis- 
cussed.) The resulting normalized numerical relations are intended to 
serve as an order-of magnitude basis for preliminary design estimates 
and comparisons. 

The response relations are derived first for a single-axis model. 
Then their applicability to three-axis design is discussed. A control 
system is postulated which decouples vehicle dynamics so that vehicle 
motions are exactly single axis. (Some advantages of such control are 
discussed in References ( 2 ) and (s).) 

The resulting control-wheel motions may be complicated by gyroscopic 
coupling due to the spinning wheels. In control to a rotating reference 
extra power is consumed also because the spin momentum of the roll and 
yaw wheels must be passed back and forth from one to the other. 

Control systems which merely damp the natural motions of stable, 
local -vertical satellites can be smaller and simpler and use less power 
but, of coiirse, furnish less precise control, (They are generally use- 
ful only in nearly-circular orbits.) 


*Based on work done for Systems Corporation of America, under con- 
tract with the U.S, Air Force, WADD Flight Control Laboratory, AF 33(6l6)- 
66'jk, and on research supported by a grant to Stanford University by the 
NASA. Part of the results in this paper are reported in Reference (l). 

** Assoc. Prof, of Aeronautical and Electrical Engineering, Stanford 
University 
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Fig. 1. 


Aerodynamic Surface. 

Aerodynamic Moment arm. 

Energy consumed, by control system 

Naperian base. Also, orbit eccentricity (introduction only). 
Vector angular momentum of a system. 

Initial angular momentum of reaction wheels. 

Vehicle moment of inertia. 

Reaction wheel moment of inertia. 

Unit vector along imaginary axis of s -plane. 

Control system gain (see Fig. 3) . 

Proportionality constant for attitude-dependent torques. 
External torque on vehicle. 

Laplace transform ofL. 

External torque on a system. 

Constant angular velocity of a local-vertical reference in 
a circular orbit. 

Power required by control system. 

Symbol representing the operation of differentiation: 

n = <i 
^ m; • 


Parameter used in Laplace transform. 

Torque exerted on vehicle by control device. 

Laplace transform of 

Time. 

Unit vectors along principal axes of vehicle. 
Unit vectors along reference axes. 

Unit impulse. 
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Time constant of controlled- vehicle response. 

Vehicle attitude deviation from the reference direction. 
Laplace transform of 0. 

A control system parameter (Fig. 3) • 

Phase angle. 

Reaction wheel speed. 

Initial reaction- wheel speed. 

Laplace transform of Q,. 

Angular velocity of reference frame. (When cd is constant 
it equals n. ) 

Forcing frequency. 


Subscripts 


1 . 2,3 

j 

V 

max. 

r 

d 

P.i’.y 


Component of a quantity along vehicle axis 1, 2, or 3- 
1, 2, or 3. 

Pertaining to vehicle . 

Maximum value . 

Pertaining to reference . 

Pertaining to damping system of last section. 

Pitch, roll, yaw. 
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INTRODUCTION 


Combination Systems . 

For common space vehicle attitude control systems a combination 
of both momentum expulsion devices (e.g. , jets) and moment'um storage 
devices (e.g., reaction vheels or gyros) are often used as actuators. 

The jets and momentim storage devices complement one another, the storage 
devices countering cyclical torques on the vehicle without loss, and 
the jets countering long-term sectilar torques by periodically expelling 
momentum from the storage system as the storage devices near spin satura- 
tion. 


Dynamic Considerations in Design . 

In combination systems the jets usually operate as stand-by devices. 
Typically, the momentum stored in the storage elements is monitored and 
the appropriate jets are fired when stored momentum reaches some percentage 
of total storage capacity. In such systems design requirements for the 
jet system are evident: (l) the total stored gas must represent a momen- 

tum capacity somewhat larger than the total secular impulse anticipated 
over the life of the vehicle; (2) the maximum torque capacity of the jets 
must be somewhat larger than the maximum torque anticipated on the vehicle — 
for example, the misalignment torque during on-orbit maneuver rocket fir- 
ing; and (3) the impulse of the jet system (moment times time) should be 
controllable to within a few percent of the total storage capacity of the 
momentum storage system (so that momentum can be expelled from the storage 
system with precision). In all of these requirements the dynamic behavior 
of the jet-controlled vehicle is of secondary importance compared with 
precision control of the jet thrust level and the thrust time. 

The momentum storage elements must be capable of storing, without 
saturation, the largest cyclical impulses anticipated plus the secular 
momentum change between momentum expulsions. 

But, in addition, the momentum storage section draws also the assignment 
of furnishing the specified attitude -control precision and speed of re- 
sponse. Moreover, it must consume minimum power and energy in the pro- 
cess. D^'namic behavior is therefore of p rime importance here. 

The present paper is concerned with this latter aspect of attitude 
control system design- -the dynamic behavior and attendant mechanical 
power requirements of reaction-system-controlled space vehicles. In 
the first section of the paper response relations are developed for a 
single-axis model of a space vehicle controlled by a reaction wheel. 

In the two succeeding sections three-axis behavior of a vehicle is 
studied for control using a set of reaction wheels. 

In subsequent papers the behavior cf gyi-c -cor .trolled systems 
will also be studied. 
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Typical Disturbances . 


Two forms of disturbance torque— an impulse and a sinusoid — 
have been chosen for the study of dynamic response. Sinusoidal dis- 
turbances commonly predominate for a vehicle in a circular orbit 
about a planet. For example, on a space vehicle whose reference is 
inertially fixed (such as the orbiting astronomical observatory) the 
aerodynamic torque will vary sinusoidally at orbital frequency and the 
gravity-gradient torque at twice orbital frequency. (At the OAO 
altitude gravity-gradient torque will predominate.) 

For a vehicle controlled to the local vertical aerodynamic torque 
may also be sinusoidal (at orbital frequency) if the orbit is slightly 
eccentric, as shown in Fig. 1 . 

Figure 1 also shows that if the orbit eccentricity is somewhat 
larger, aerodynamic torque on the vehicle is essentially an impulse. 

The reason is, of course, that atmospheric density falls away so 
rapidly with altitude that in an eccentric orbit the vehicle spends 
only a small portion of its orbit period in the effective atmosphere. 

For intermediate values of eccentricity the aerodynamic torque may 
be represented by a Fourier series, involving the first few multiples 
of orbit frequency, so that, again, sinusoidal response is of interest. 
(The Fourier series approach has been employed, for example, by 
Schrello in Reference (4).) Soiirces and magnitudes of disturbing torques 
have been discussed in a number of papers (c.f. References (4), ( 5 ), 

(6)). Certain of their effects will be discussed further in Reference 

(3). 


Other instances of Impulsive distiorbances may occur (e.g., meteorite 
impact). More generally, impulse response furnishes a good starting 
point for analysis of response to random disturbances. In the present 
paper sinusoidal and impulse dist\irbances of typical size will be used 
to illustrate the form and order of magnitude of salient response rela- 
tions . 


SINGLE AXIS CONTROL 


Figure 2 offers a single -axis model of a space vehicle with a 
reaction wheel for attitude control. It is shown in the next section 
that this model can represent the actual situation quite well in many 
situations of interest; namely, when the attitude stays near a reference 
which is inertially non-rotating, wheel speeds are slow, wheels are 
accurately aligned to vehicle principal axes, and the vehicle is 


* 

Calculations for Fig. 1 are by Roger Bourke, Stanford University. 

** 

Parts of this section are repeated in Reference (2). 

The variables are attitude deviation, 6, and wheel speed, I — 

is an external disturbing torque, and | ^ is the internal control torque 
between the wheel and the vehicle. 
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aerodynamic torque on a hypothetical local level satellite 


syiimetrical with respect to external torques. In such situations dynamic 
behavior about each axis is decoupled and may be treated independently. 
Moreover, as will be shown, single-axis results furnish a most useful 
reference for three-axis studies. 

In the present simple study it is ass\mied also that control will 
be "stiff" compared with physical restoring torques, such as gravity- 
gradient or aerodynamic 0-dependent torques, and these torques will be 
omitted at this point. Their effect is discussed at the end of the 
present section, and also in the next section. 


* 

fiynamic Equations 

Three equations of motion can be written for the system of Fig. 
2, any two of which constitute an independent set in the variables 0 
and fi. Newton's Second Law is applied, in turn, to the systems of 
Fig. 2(a) (vehicle plus wheel), 2(b) (vehicle only), and 2(c) (wheel 
only) . 

From Fig. 2(a), for the vehicle plus wheel: 


(I + J)0 + JS2 = L 


or, after Laplace transformation: 


(I + J)s^0 + Jsil= L (la) 

From Fig. l(b), for the vehicle only: 

Is^e = L + T (lb) 

c 

From Fig. l(c), for the wheel only: 

Js(iX+ se) = -T (ic) 

c 


* 

The equations of motion for a much more general situation have 
been derived by Roberson, cf.. Reference (t) 
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FIG. 2 SINGLE-AXIS MODEL OF VEHICLE WITH REACTION WHEEL 



In Eq. (l) T is the control torque, applied by the control sys- 
tem between the wheel and the vehicle. SI is the angular velocity of 
the wheel relative to the vehicle. Note that friction torque and 
back e.m.f. between wheel and vehicle are assumed zero in this study. 


Control System . 

A simple control system is shown in Fig. 3 for controlling the 
attitude of the vehicle in Fig. 2. An attitude sensor (e.g. , a star- 
sight system) is assiamed to report 9 . 


The upper part of the block diagram in Fig. 3 is taken directly from 
Eq. (lb) . To this is added a simple control loop employing proportional 
plus derivative control. For simplicity the control gain is set so that 
the system characteristic has two equal real poles, as given by the root 
locus picture in Fig. 3b. The gain required for this is: 



T 


in which t is the final time constant of the system. 


With this gain setting, the overall system transfer function is 


2 

® = I_ 1 
L I (ts + l)^ 


( 2 ) 


Response to Initial Attitude Error . 

The response of the system of Fig. 3 to an initial value of 9 is: 

9 = e(o)(^l + ^ (3) 

Response to initial 6 is plotted to a linear scale in Fig. 4a, 
and to a semi-log scale in Fig. 4b. Choice of t for adequately fast 
recovery from initial errors can be made from Fig. 4. For example, if 
the initial misalignment must be reduced by a factor of 100 in I 5 OO 
sec. then, from Fig. 4b, 6.6t must equal 1500, and the required system 
time constant is t < 250 sec. 


Response Relations for Impulse Disturbance , 
(a. ) Attitude Excursion 


If U is an impulse of magnitude i = L At, as shown in Fig. 5a^ 
then its Laplace transform is L(s) = f, and €he system response is given 

t>y 
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(a) System Arrangement 


FIG. 3 ATTITUDE CONTROL SYSTEM 
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iT t 
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This relation is plotted in Fig. 5^3. 

The maximum deviation of the vehicle occurs at t = and has the 
value : 


e 

max. 


£t 1 
I e 


(5) 


For purposes of numerical calculation, for specific design problems, 
it mi^t be convenient to have available a set of response curves for 
various values of t. Such a set is given in Fig. 6, from which not 
only can the maximum vehicle deviation (from its reference attitude) be 
read, but also the length of time the deviation will exceed any specified 
value . 


For example, suppose the appropriate control-system time constant is 


whose peak 
is 


2^*^ gm cm'- 

attitude excursion, under an impulsive disturbance of ^ x dn cm sec, is 
to be held under 10 mr, and- -further — whose misalignment may exceed Imr 

for no longer than 200 .sec after the same distiurbance. The first specif ica- 

°max iNr^ /- 

tion gives ■ ^ t-tj-t:- = 50, for which, from Fig. 6, t must be 


to be chosen for a vehicle of moment of inertia 


made less than about l40 sec. But the second specification requires that 
be larger than = 5 for no longer than 200 sec which, 

from Fig. 6, means t must be less than about 50 sec. 


i/I 


(b.) Control Torque Required . 

From Fig. 3 and Eq. ( 2 ) wheel control torque is to be related to 
vehicle disturbing torque by the transfer function: 

T = - ^ (2ts + 1)0 = - L (6) 

(ts + if 

from which, for L an impulse of magnitude £ , the time response is: 

t 

T^(t) . - ^ (1 - |-) e' " (7) 

This response is plotted in Figure J. Note that the control torque has 
its peak magnitude at t = 0, the instant at which the disturbing impulse 
occurs. Its value is: 


The symbol is used to abbreviate "xio’^." For example, 
means 2x107. 
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time ,t ^ 

(a) Impulsive Disturbing torque 



(b) Vehicle Response 


FIG. 5- VEHICLE RESPONSE TO IMPULSIVE DISTURBANCE, 
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FIG. 6 DESIGN RELATIONS 
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(Ta) 


A. more direct way to derive the desired control is to solve Eq. (la) 
for XL or (lb) for after substituting for L from the desired- 

response function, ( 2 ). The result is : 


T = - JsXl= + ls^0 - ^ (ts+1)‘^0= - %(2ts+1)0 , 

C d d. 

T T 


XI = 


I 2ts+1 
Jt^ s 


T^= -~2 (2ts + 1) 0 , 

T 


the last expression being the same as (6). 


(c.) Wheel Motion Required . 

From Eqs. (ic) and (6) or ( 2 ) and (8a)the relation between the 
vehicle distiurbing torque and control wheel motion is given by: 




[2ts + Ij 


Js (ts + 1) 

Again, if the disturbance is an impulse, the time response is: 


J2(t) =i|^l + ( L-i) e , (9) 

which is plotted in Fig. 8 for initial speed, = 0. The wheel, of 
course, rotates in a positive direction (by the convention of Fig. 2), 
so that it, rather than the vehicle, will have the momentum change required 
by the application of torque impulse £ to the overall system. 

(d.) Power and Energy Required . 

The minimiom value of power required by the control system will be 
that needed to accelerate the reaction wheel, assuming no f ri ction, and 
assming that no power is dissipated in electrical losses . This 

power is simply the product of control torque and wheel speed relative 
to the vehicle. (This is a dot product — that is, it involves the sign 
of both quantities). 


■^2 J 

Term Js 0 is dropped in (la) on the groiinds that — < < 1. 

** 

Only this limiting value of mechanical power is considered in the 
present paper. 
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FIG. 7 


m 


FIG. 8 



CONTROL TORQUE RESPONSE TO IMPULSIVE DISTURBANCE 



WHEEL SPEED RESPONSE TO IMPULSIVE DISTURBANCE 
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Thus, for an impulsive disturhance , power is given by the product 
of -Te(t) fromEq^.(7) and 12(t) fromEq.(9)- (The minus sign on 
is necessary because Eq..(7) gives torque on the vehicle while the present 
calculation is concerned with torque on the wheel. Note that the sign 
convention in Fig. 2 considers torque on the wheel and wheel velocity 
positive in opposite directions.) For the special case that = 0 
the power is : 


2 / 

p(t) = 


T L 




- 2 ^ 

T 


(lOa) 


P(t) is plotted in Fig. Maximum power is required at about t = — 

and has the value : 

P = .66 ^ 
max tJ 



Note that the power, as plotted, actually becomes negative after 
time t = 2t, when the wheel is being decelerated as it approaches its 
final speed. This result is based on the premise of recovering energy 
from the spinning wheels when they are to be slowed down, and storing 
the energy (for example, in another wheel spinning in the opposite 
direction). The feasibility of such energy conservation is discussed in 
Appendix B. If energy cannot be recovered then, of course, the magnitude 
“only^ should be read in Fig. 9* 

The total energy consumed is foum'by integrating Eq.(lOa) or find- 
ing the area under the curve in Fig. ^a. The value obtained is: 

.2 


E = .5 


max 


(11a) 


assuming, again, that negative values in Fig. 9 represent recoverable 
power. Energy consumed would, of course, be slightly larger if no 
recovery were assumed. 


If a wheel is spinning with some high initial speed, K , when 
the vehicle is disturbed (e.g., by impulse, i), then a much higher power 
level is required to accelerate the wheel. In this case the power is 
given by an amended version of (lO) (obtained by multiplying ( 7 ) by ( 9 )): 


P = 


2i 

Jt 


(1 - 2 t ^ 


t 

T 




-2 


2 ^ , - - 
+ -^(1 ^ ^ (10b) 


in which the last term may be much larger. When this is so the power 
consumed will look like Fig. 9h, and the energy will also be much larger: 

E = (lib) 

That is, the torquer must expend much additional energy to chase the 
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(a) Zero Lfiitial Wheel Speed 



(b) High Initial Wheel Speed, 


FIG. 9 POWER CONSUMED IN RESPONDING TO IMPULSE 
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spinning sheel. This energy is larger than the zero-initial -vheel-speed 
energy by the ratio 


ia JC h 



in which h is initial wheel momentum. 


Response Relations for Sinusoidal Disturbance . 


Suppose that the system just described is dist\irbed by a sinusoidal 
external torque: 

l_ = L cos w t ( 12 ) 

o f 

Then, from ( 2 ) (using, for example, the steady-state portion of the 
Laplace transform) : 


max I y 


in which 


L 

““O 


max 


(ts + 1) J 


s=ju 


( 13 ) 


and 



L(ts +1) J 




Usiially u will be orbital frequency, while t will be shorter than 
i/ 27 T times the orbital period, so that < 1, and the relations are 

simpler: o 

I 

(li^) 


L/ 

o 


For tu)^ <1; ® ~ I 


cos <*>^t 


The corresponding values of control torque and wheel speed are, 
from (6) and (8 ) \ ^ 


T = - L cos w t 
c o f 


For TU)^ < 1 ^ 


k . 


£2 = sin w t 


( 15 ) 


(16) 


Then the power required is: 
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L ^ sin 2w t 

which is plotted in Fig. 10. 

It is important to note that the power has no secular component, 
which might have resulted from rectification when two sine waves were 
multiplied together. There will never he rectification (for the ideal, 
no-damping model assumed here) because SI is always the exact integral 
of : they will always be 90° out of phase . 

If energy can be recovered as the wheels are driven, then the net 
energy per cycle will be zero. But if, on the other hand, a drive -then- 
brake method is used, then the energy consimed will be the shaded area 
in Fig. 10 (braking involves no power) : 

L ^ . 

E = per cycle = 2 *^/°-'x(l8) 

jco 2 f 


As an example, suppose a vehicle having moment of inertia 
I = ^0 gm cm^ is controlled to an inertial reference with wheels of Inertia 
J = 2 \ gm cm^ and with a time constant of 100 sec. Suppose, the vehicle is in 
a slightly eccentric orbit such that the sum of aerodynamic and gravity- gradleni,, 
torque is sinusoidal, at orbit frequency, with peak magnitude 1 v dn cm. 

Then, from (l4), (15), (16), and (it)^ the peak values of vehicle attitude 
excursion, control torque, wheel speed, and power required will be, 
respectively. 


max 


Lt" 

o 


(i\^)(iV) 


,005 rad. 


T = L = 1^ dn cm 

o 

L l\^ 


c o 

max 


’ (At (iV^) 


= 500 rad/sec. 


P = SI 

max 2 c max 2 
max 


_ 1 (i\**) (;V) _ 




= 25 watts 


(The factor l'^ converts dn cm/sec. to watts.) 

Without energy recovery, the energy used per cycle, by (18), is: 

E = = — : — r — = 500 watt-sec. 


(2\^) (iV^)^ 
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FIG. 10 POWER CONSUMED DUE TO SINUSOIDAL DISTURBANCE 
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In cases where too is larger than 1 (i.e. , the system time con- 
stant can he made longer than an orbital period j or a high-frequency dis- 
turbance is involved) the above quantities may be attenuated^ as indicated 
by the functions of s in ( 2 ), (5) and (8), This is discussed further 
in the last section of this paper. 


Effect of Physical Restoring Torques . 

Sometimes, if the system time constant is to be very slow, physical 
0 -dependent torques --particularly gravity-gradient or aerodynamic torques — 
may be important factors in the dynamic characteristices of the system. 

In such cases the external torque on the vehicle consists not only 
of the independent component, L , in Eq. (la), but also of a 0-dependent 
component. In the case of a local -vertical satellite, and for the very 
small 0's we are considering, this torque may be written as a linear 
function, -k0. Then Eq. (la) becomes 

Is^e + JskQ= L - k6 


or 



+ k)0 + Js^l.= L 


(19) 


In applications where positive, precise control is indicated, the 
presence of such additional torques may be taken care of by having the 
control system cancel them out. That is, the term k0 in (19) would 
be countered by amending the control torque equation, (8a), as follows: 


Jsil,= 


^ ( 2 ts + 1) -k 


0 


(20) 


With this control the attitude response of the vehicle will be exactly 
as given by (2), as can be seen by substitution of (20) into (l9)» 

Control torque and wheel speed would, of course, be altered by the addi- 
tional control term k0. As a practical matter, however, the magnitudes 
of physical torques will be small, compared to control torques involved 
in precise attitude control, so that only small percentage changes in 

y and would be involved, 
c 

Occasionally only light attitude stability is required and only 
approximate attitude angle must be maintained. In such cases it may be 
possible to use physical restoring torques to supply the main stabilizing 
torque, with the attitude control system functioning only to supply 
damping. This situation is depicted in Fig. 11a, which is similar to 
Fig. 3, except that the vehicle dynamics now include the physical re- 
storing torque (e.g. , due to gravity gradient). The corresponding root 
locus picture for the case that k is positive (stable restoring, torque) , 
is shown in Fig. 11b. Without the control system, the vehicle would 
oscillate indefinitely at its own natirral frequency. The control system 
supplies damping, so that natural motions of the vehicle are damped out 
in time. 
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( b ) k Stabilizing ( c ) k Destabilizing 


FIG. 1 1 CONTROL WITH PHYSICAL RESTORING TOROUE 
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The use of passive stabilisation (with control system damping) is 
probably confined to those applications in which the orbit can be 
counted upon to be nearly circular (e.g. , eccentricity < .OOl). Otherwise^ 
as some studies have shown, (c.f. References (8) and (9)), the gravity 
gradient torque may actually cause dynamic instability for some initial 
conditions of the vehicle. 


THREE -AXIS CONTROL 


I^rnamic Equations . 

The equations of motion for a space vehicle with internal moving 
parts have been written elsewhere for quite general circumstances. 

(c.f. Reference (t).) However, it will be convenient to rederive the 
dynamic equations briefly for the special case of interest here. 

_ _ Consider the space vehicle shown in Fig. 12 with principal axes 

1, 2, and 3 fixed in the vehicle, and assume the special case that the 
three control wheels have their spin axes aligned exactly along the 
principal axes of the vehicle, as shown. The attitude of the vehicle 
is to be controlled to the reference axes Ij., 2^., and 3^., which are 
also an orthogonal set. (in general, the reference axes may be 
inertially fixed or may rotate in some specified manner). Let the 
orientation of the vehicle with respect to the reference axes be defined 
by angles 6 , 0 , 0_, as shown, (it is assumed that 0's are always 
small. ) 

Next, assume further that the reference system is either an inertial 
one, or is orbit -oriented to the local vertical and the orb^t axis. 

In the latter case, to simplify the equations, we let the 2p axis co- 
incide with the orbit axis and the 5 ^ axis coincide with the vertical, 
£0 that the reference system has an angular velocity, say about the 
2j, axis only. (Except for a perfectly circular orbit, w will vary 
with time.) The equations written for this reference system can be 
quickly specialized to the inertial reference case by letting w equal 
zero. 


The system of Fig. 12 has three degrees of freedom to locate its 
mass center, plus six more to define motion with respect to the mass 
center (three for the rigid vehicle, plus one each for the wheels 
relative to the vehicle). It can be shown that, for situations of 
interest, motions about the mass center will have a negligible effect on 
motion of the mass center. For present purposes, therefore, we can con- 
sider variables defining motion of the mass center as independent, and 
write Just six dynamic equations for motion about the mass center. We 
shall write them using the three vehicle attitude angles, 0, plus the 
three wheel speeds, Q as the dependent variables. (in each case we 
write the well known extension of Newton's Second Law for motion of a 
system of rigid bodies about their mass center.) 
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For the first set of three equations we elect to write the law for 
motion of the entire system — vehicle plus three wheels: 


I system ^ veh. 




= M 


Details 

of 

the derivation are given 

in Appendix A. 

The result is: 

+ 
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0 
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+ 

J 

0 

P 

0 

n = M - 
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-U) 

0 

P__ 


0 





in which p represents the operation — , the components of 9 

Q.'C 

are attitude angles 9-^, 02 ^ f -he components of K are the 

angular velocities of the three wheels, ^2^ hhe components of 

M are the hody^axis components of external torque, M^, M , M^. (j has 
been lumped with I in the 9 terms . ) 

Often external torques are attitude dependent, as discussed in the 
preceding section. Typically, aerodynamic torque will have a specific 
time variation if the vehicle remains always oriented exactly to its 
reference axes; but if the vehicle deviates, then the torque may vary 
also as a function of attitude kngle 9, In the above equations, there- 
fore, it will be convenient to divide total external torque M into a 
component, L > which is independent of attitude plus the 0 -dependent 
components , which we assume can be written as linear functions of 0 
because we shall confine ourselves to very small 0. Then the above 
equations become: 



JBjP + 


-J2^p - Ja^cD + 

+ '"22 

Jft^p - 
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Note that the h’s can he either positive or negative, and can he time 
varying. Fiirther, although the k’s due to gravity gradient will appear 
only on the diagonal of the matrix, because 1, 2, 3 are principal axes 
in Fig. 12, other torques — for example the aerodynamic k's — may not he 
symmetrical (e.g. , because of paddles unsyrametrically deployed): hence 
the off-diagonal k’s in Eq. (2l) . 

In this paper we wish to study some response relations for the 
system under the simplifying assumptions that the equations are linear 
with constant coefficients. In some cases these assumptions will he 
quite accurate. For examgjle, if changes in the wheel speeds are to he 
small in a given respopse, then the perturbation technique can he used 
to linearize the (J2)& terms. If, further, the orbit is circular, 
then w is constant, (if the reference is an inertial one, w is 
zero.) Under these conditions Eqs. (2l) become linear with constant 
coefficients. 

In other cases the assumption of linearity and/ or constant coeffi- 
cients can serve only as a first approximation, the utility of which 
must he evaliiated in each specific case. 


The second set of three dynamic equations is obtained by writing 
Newton’s law for each of the wheels about its spin axis; 

j(a, + e.) = - T 

J 0 =j 

iri which j = 1, 2, or 3. (As before, the minus sign occurs because 
^ , /I torque between the wheel and the vehicle, is considered positive on 

the vehicle and negative on the wheel.) In most cases of interest, 9 
will always be very small (of order j/t) compared to , and will there- 
for be dropped: 



( 22 ) 


Control to an Inertial Reference . 

Note, first, that Eqs.(2l) can be linearized and completely decoupled 
on the basis of a specific set of assumptions. These are (l) that the 
reference frame is an Inertial one, (2) that the wheel speeds are 
initially near zero and do not change much during the motions of interest, 
and ( 3 ) that the vehicle has aerodynamic, as well as inertial, symmetry 
with respect to its principal axes so that cross -torques, kj^j, are all 
zero. These assumptions produce, respectively, the following simplifica- 
tions in ( 21 ): ( 1 ) 03 = 0 , ( 2 ) (Ja)e = 0, and ( 3 ) k^. = 0. Then all off- 

diagonal terms disappear, all on-diagonal terms are linearized, and (2l) 
may be Laplace transformed (if desired) giving: 
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+ Jsfl= L (23) 


V "•'33- 


Each of Eqs.(23) is an independent linear equation analogous to 
Eq. (la) . Correspondingly, Eq. (22) is analogous to (ic). Therefore, 
all of the results of the section on single axis control apply verbatim 
to the system described in the preceding paragraph. 


Consider next the sitmtion in which we remove special assumptions 
( 2 ) and ( 3 ) above, but retain assumption number (l) (that the reference 
is an inertial one). In this case, since we^are permitting the wheels 
to have high initial speeds, the terms (Jfi)0 in (2l) will be of first- 
order importance. However, if the changes in wheel speed during a 
period of interest are small compared to the initial speed, then the 
equations can still be linearized by using the perturbation technique. 
That is, we assume that each wheel speed consists of a constant initial 
valueplus a small perturbation: 


fl = + C. 

j JO J 


Then, after Laplace transformation. 


Eqs.(2l) become 
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in which h = JQ . The gyroscopic coupling terms are very much in 
evidence . 

It is shown in Reference ( 2 ) that there is considerable advantage 
in choosing the control equations for a coupled system on the basis of 
desired performance and, further, of specifying that the responses about 
the three axes be decoupled from one another. Specifically, let it be 
required that a disturbance about the principal axis of the vehicle pro- 
duce a critically-damped response about that axis (as in the single- 
axis case) and no response about either of the other axes. That is, 
the response is to be: 
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which has already been plotted in Fig. 5 for L an impulse. The rela- 
tion between r and ^max is given ^n more detail in Fig. 6. Equation 
(24) can now be solved directly for , substituting for L from (25)'- 
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Thus, while the vehicle response to a disturbance is Identical with 
the single-axis case, the wheel-speed response involved is now more 
complicated. It can be obtained directly in terms of disturbance L by 
substituting (25) into (26) and solving for\Qy- The result is: 
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(27) 


Each of the terms in (27) is of the form: 


C L(s) 


If L is an impulse, then the time response corresponding to each term 
in (27) will resemble, in both form and magnitude, the single-axis re- 
sponse given by Eq. (9) and plotted in Fig. 8, as is evident from a compari- 
son of terms in (27) with (9)* 


The point indicated by Eqs.(27) is that with high Initial wheel speed 
the wheels will be two or three times as busy and will respond to inputs 
about all three axes. 

, 2 

The terms in (27) represent the relative "stiffness" of the con- 

trol system and the natirral restoring (or destabilizing) torques. In the 
usual case, the control system will dominate, so that these terms will be 
small. We therefore drop them at this point. (Their inclusion adds a little 
complication, but is straightforward, as was discussed in the Single-Axis 
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Sectiono The special case of a control system which damps natural motions 
of a local level system is discussed in the last section of the paper.) 

h.T 

Terms represent the coupling strength, of course, and may he as 

large as 1 for high-capacity wheels. 

The control torques involved in accelerating the wheels as specified 
in (27) are available from (22): 

T = - Jsfl 
c 


in which the components of a are given by (27). The form of torque 
response will therefore be like that given by Eq, (7) and plotted in 
Fig. 7. 

To obtain power consumption we recall that we have selected the case 
in which change in wheel speed is much smaller than the initial value, 

The power will therefore be given, closely, by: 


P = - T * fir 


That is; 
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(28) 


Two forms of time response are obtained from (28), namely: 


(2 -^) 


, and (1 - — ) e ^ 


These are plotted in Fig. 13 over time periods of interest. 

The total energy consumed by the control system during its response 

to a set of impulses, i , i , i , about the three vehicle axes, is given 
by 12 3 


E 


f>CO 

^t =0 


[P^(t) + p^Ct) + P^(t)] dt 
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From the resiilts of Fig. 15, the total energy can be written explicitly for 
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( 2 ) energy is recovered (areas in Fig. 13 are added with signs as shown): 

E . ! (29b) 


The last result suggests the striking possibility that energy con- 
sumption may be reduced to the conservative, single-axis level (Eq. (lib)) 
without storage, by effecting the appropriate transfer of energy from one 
control wheel directly into another, the proper transfer program to be 
specified by the computer. 

In the absence of energy recovery the total energy, (29a) will be 
somewhat higher than for the single-axis case. 


Response to Sinusoidal Disturbance . 

If the time constant of the control system is short compared to the 
period of a sinusoidal disturbance to which the system is subjected, 
then — as discussed under single-axis response — the response of the system 
is greatly simplified because terms can be dropped compared to 1 

in expressions like ( 25 ) • Commonly the disturbing frequency is orbital 
or twice orbital frequency and the system time constant is many times 
faster, so that this is a good approximation. 

If the distiurbance is 

L =r L cos “ t , 

o I 

and if < 1, then the attitude response (which is single-axis be- 

cause of the decoupling control) is, from ( 23 ) • 
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Lt2 

6 . = — cos w t 

about each axis (which is the same response as (l 4 )). 

The wheel speed and control torque- -and hence power and energy- - 
also revert to the single-axis form because all coupling terms in (26) 
contain . Thus SI, T^, F, and E are given by single-axis rela- 

tions (16), (15) ^ (T?)^ and (18), respectively. 




COHTROL TO A ROTATING REFERENCE 


If the reference frame is rotating, at angular velocit;^- then 
two conditions are necessary to linearize Eqs. ( 2 l) : (l) 6 must be 
small compared to (2) w must be small compared, for example, to 

Y If both these conditions are met then (21) can be linearized and 

Laplace transformed, giving (with t*> = n) : 
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( 30 ) 


provided only that the vehicle has symmetry with respect to external 
torques. (if this is not true, then the 6 matrix will have k. ,*s 
in place of its O's. and the analysis will be somewhat more compiicated, 
but straightforward.) ■ , 

In studying the response to sinusoidal disturbance at orbital 
frequency condition (l), above, reduces to 

n 6 « n 
max 

which is met by the restriction to small 6 . (0 is deviation from the 

reference axis system.) For faster transient motions condition (l) may 
be met only approximately, and must be assessed in each case. 


Except for terms kj^j (e.g. , aerodynamic coupling due to non- 
symmetry) • In the unusual case that these are large enough to cause 
serious coupling the present analysis must be modified accordingly. 
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Condition (2) is met exactly, of course, only by a circular orbit. 
Deviation from condition (2) depends on orbit eccentricity, and must also 
be evaluated in each case to determine the applicability of linear analysis 
based on (30) . (if deviation is too large then, of course, linear analysis 
can give only a rough approximation to the behavior, and non-linear tech- 
niques or computer solutions are indicated. ) 

It is convenient to divide systems for controlling a local-level 
vehicle into tvo classes: precise control systems, and systems for damping 
natiiral motions. 


Precise Control System. 


If precise, relatively fast control of the vehicle represented by 
(30) is required then (l) the torques involving n and k in the 0 
matrix of (30) will t\rrn out to be small compared with the control 
torques, and will have negligible effect on the dynamics, but (2) the 
n terms in the Q. matrix have an important effect involving initial 
wheel speed. 


To show these two points most directly we choose a decoupling con- 
trol system, such that vehicle motion will be given, again, by: 
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in which pitch and roll -yaw have heen decoupled, and initial wheel speeds 
have heen included. 

Fo r precise control t will he very short compared to either 
or l/n ^ and corresponding terms may be dropped from (3l), 
proving point (l). Then (Sl) becomes: 
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(32) 


From Eqs.(32) the pitch wheel motion in orbit behaves just as in 
the single-axis case (compare with (8)), but the roll and yaw wheels 
have now also a sinusoidal motion at orbital frequency. This motion will 
persist — if there is an initial wheel speed- - even in the absence of 
disturbances . This is point ( 2 ). 


Physically the reason is that the momentum vector stored in the 
wheels must be maintained fixed in space and, therefore, must be passed 
from one wheel to the other as the vehicle rotates at angular velocity 
n. This requires control torque, of coiorse, and therefore power. 

The power required is obtained by combining (22) with (32). For 
no disturbance and C 2 (°) “ example the amount is 


P = Jn C^(o)^ 


sin 2nt 
2 


( 33 ) 
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which may represent a substantial proportion of the total mechanical 
power required. 

The above situation is investigated in greater detail in Reference 
(3). In particular it is shown there that if decoupling control is not 
used the momentum exchange must be accompanied by an attitude error which 
may be very large. 

It should be pointed out that the only way to avoid the requirement 
to accelerate and decellerate in a local-level satellite is to use spin- 
ning members which are so mounted that their spin axes are free to remain 
inert ially fixed- -as with the reaction sphere. 


System for Damping Natural Motions . 

From (30) the possibility is evident of using physical restoring 
torques (e.g., k's) to stabilize a local -vertical satellite. The most 
commonly used torque is the gravity gradient. 

The natural dynamic characteristics of a rigid vehicle in a 
gravity-gradient field have been described carefully elsewhere (c.f. 
References (8) (9) and (ll).) To summarize briefly, the corresponding 
values of the k's in (30) are, for a symmetrical vehicle: 


^11 ^^ 2 "^ 3 ^ ^ 


(roll) 


^22 ^ 
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(pitch) 


^ 33 = ° 


(yaw) 


in which Kg is a gravitational constant. Further, for a circiilar 
orbit, K /r^ = n^ , so that (30) becomes 

s 








I3S +(I2-Ii)n 63 
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in which is the inertia of the damping wheel. 
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Evidently, with the proper relative magnitudes of the I's the 
vehicle will he neutrally stable due to gravity-gradient alone --i.e., 
with the wheel speeds zero. That is, it will perform undamped oscilla- 
tions about all three axes at simple multiples of orbit frequency. 

For example the celebrated dumbbell configuration oriented so that 
= 0, will oscillate at frequency v3 n in pitch and 2n in 
roll. "^More generally, requirements on the I’s are that either 
I > I > I or part of the region I > > Ip for stability. Dynamic 

c^racteristics of rigid satellite vehicles are aescribed comprehensively 
in Reference (ll) . Note that yaw stability is achieved because of roll- 
yaw inertial coupling produced by the rotating reference frame. Basically, 
this is a form of gyrocompassing. 

From ( 34 ) it is easy to see how light damping of a naturally- stable 
vehicle can be accomplished. In pitch, for example, we can convert (34) 
to 

I^Cs^ + 2:;^ + u,2) ©2 = ( 35 ) 

d ^ 


by letting 






(36) 


(‘*’ is an abbreviation for the vehicle natural frequency in pitch, 
P 

1 

13 2 

y — 3n , T, is the time constant of the damping envelope, and 
2 

is damping wheel inertia), 
d 
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From ( 35 ) the response to an impulse disturbance, L-g = , will 
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Corresponding wheel speed will be 


2i 


= 


2 J T 

d d pd 


sin ^ ,t 
pd 


(38) 


and control torque: 
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for ^ < 1. 
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(C is the damping ratio ^ .) 

d p 

Power is given by the product of (38) and (39): 
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Comparison of (37)^ with shows that for the li^t-damping sys- 

tem attitude excursions will be larger than for precision control by 

iZi 


Comparison of (38) with (8) 
speeds, the damping wheel may be 


shows that, for the same control -wheel 
made smaller by: 
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d pd 
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In this case comparison of (39) and (Uo) with ( 6 ) and (9) shows that 
control torque and power will both be lower, by x ,/ t, for the damping 
system: 
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(41) 


Equation (ii-l) is a convenient relation for preliminary design studies. 


In roll-yaw the decoupling principle is again useful, and we can 
specify response by 
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in which w _ —§—-5. u ^ n^ • (The same 

’ y I3 

is chosen in all three axes in this example.) 

Then the wheel program required will be, from (34): 


(42) 
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From (42) the attitude response in roll and yaw will he the same as 
in pitch, (37). 

The wheel motion (and therefore, also, torque and power) in roll- 
yaw will he somewhat more involved than in pitch because of the off- 
diagonal terms in (43). Substituting for 0 from (42) the response 
of wheel 1 to an impulsive torque (about an axis in the roll-yaw plane) 
is given by: 
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The expression for is similar. The first denominator quadratic 

leads to a time response similar to (38). But, in addition, the second 
denominator quadratic leads to a component of wheel speed which is 
sinusoidal at orbital frequency, just as there was in (32). Again, 
this is necessary to keep the wheel momentum vector fixed in space when 
there is no disturbance, and cannot be avoided in control to a local- 
vertical reference with the system of Fig. 12. (Note that when n=0, 
the above equation becomes just like (36).) Moreover, for the present 
system in which and Wy are slow, of order n, the peak power re- 

quired to accelerate and decelerate the wheels during the orbits follow- 
ing a single impulse will be as large as that required to damp the transient 
attitude motions produced by the Impulse. 


Sinusoidal Response . 

For a damper system the sinusoidal response relations are more in- 
volved than for a tighter control system simply because the forcing 
frequency is no longer much lower than the natural frequency of the system, 
but, on the contrary, may be somewhere between ^ the natural frequency 
(for disturbance of a dumbbell in roll at orbital frequency) to several 
times the natural frequency (e.g. , for disturbance of a more nearly 
Iso-lnertial vehicle.) 

For a gravity-stabilized vehi cle it is a ssumed that the natural 
frequency of the vehicle — e.g., t ^ in pitch — will be at least 

ys n 
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a sizeable fraction of orbital frequency. Otherwise gravity stabiliza- 
tion would be too weak to be useful in the presence of other torques. 

To illustrate the relations, the frequency response of a vehicle 
in pitch is shown in Fig. l4. Accompanying plots show wheel speed and 
torque response. It is seen to be advantageous to avoid resonance, if 
possible. 


Eccentric Orbit . 

It must be emphasized again that the present analysis assumes a 
near-circular orbit. When the orbit is eccentric the gravity gradient 
may actually be destabilizing, so that more positive control is required 
and the technique of merely damping natural motions cannot be used. 




(b) Control Torque 


FIG. 14 RESPONSE OF DAMPING SYSTEM TO SINUSOIDAL DISTURBANCE IN PITCH. 




SUMMARY 


In the common attitude control system which utilizes a comhination 
of momentum-expelling and momentum- storage mechanisms the latter draws 
the job of furnishing control precision and speed of response with 
minimum power. 

Design of a reaction wheel system for required response is aided by 
basic response relations for initial error and impulsive and sinusoidal 
disturbances, as illustrated in Figs. 4 through 10 and Eqs. (l3^) through 
(l8) of the present paper. Specifically, controlled-systera time' constant 
can be chosen from allowable attitude escursion in Fig. 6, for an impulsive, 
or Eq. (l4) for a sinusoidal disturbance; or from allowable Initial re- 
covery time in Fig. 4. Torquer capacity and wheel size can be determined 
from Figs. 7 and 8, for an impulsive disturbance, ^or from Eqs. ( 15 ) and 
( 16 ) for a sinusoidal one. Peak power and energy are g^iven by Fig. 9 
and Eq. (ll) for an impulse, or by Eq. ( 17 ) and ( 18 ) for a sinusoidal 
disturbance . 

While the above relations are derived for the single axis case, 
their applicability to three-axis design is established. For precision 
control to an inertial reference the postulation of decoupling control — 
which offers several advantages — results in vehicle response which is 
identical with the single -axis case, but control wheel activity may be 
greater by a factor of 2 or 3, due to gyroscopic coupling, if the initial 
wheel speeds are high. Even so, if a technique can be devised for re- 
covering mechanical energy, then energy consumption may be reduced to 
the single -axis level. The possibility of mechanical energy conservation 
is discussed in Appendix B. 

For precision control to a rotating reference the single-axis vehicle 
response may still be achieved, but control activity is now complicated 
by the continuous angular velocity of the vehicle. In particular, spin 
momentiim must be transferred back and forth between the roll and yaw 
wheels--at significant expense in power — to keep their momentum vector 
fixed in space. An accompanying attitude error can be avoided only by 
using decoupling control. 

A control system which merely damps the natural motions of a local- 
vertical satellite can be used if precision control is not required and 
if the satellite has a stable configuration and is in a nearly-circular 
orbit. Then the response of the vehicle will be larger and slower, roll 
and yaw motions being generally coupled, but smaller wheels can be used 
(Eq.(39)), and power consiomption will be lower by the ratio of damping 
time, in the damper system, to the control time constant of the precision 
system. 


* 

Only mechanical power, required by an idealized torquer, is 
calculated. 
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APPENDIX A: EQUATIONS OP MOTION FOR SYSTEM OF FIG. 12 


Refer to the assumptions stated under "THREE- AXIS CONTROL". 

The angular velocity of the vehicle-fixed coordinate system 1,2^3#' 


— veh — veh/ref —ref 

n = 0. ' + a 


(A-1) 


For the assumptions given, that the 0's are small, and that the refer- 
ence system rotates about 2^ only, (A-l) can be written: 

—veh • 77 

n a; 0 + 2 (D 
r 

« 1(0 + oi 0 ^) + 2(02 + od) + 


(A-2) 


That is, products of Euler angles and their derivatives are neglected. 


The angular momentum of the vehicle is given by: 


= [I] 


i^ 0 0 

0 I2 0 


—veh 
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3 


and the angular momentum of the wheels by: 


Z = J + a] 


in which 2 is defined (as in the text) as a vector made up of the wheel 
spin velocities: 



Then, 

T” 

i .e . , 


with + Z H^, we can write: 

^syst ^ fsyBt ^ -veh ^ jByst ^ 

30 — _ 

H means the derivative of.H a_s se^ by an 
the derivative of H with T = 2 = 3 = 0 . 
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* 

M (A-3) 

observer on the vehicle 



with respect to vehicle coordinate system 1,2,3* again, one drops 

the small products of Euler angles and their rates, and also compared 
to 1, the resiilt is : 
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FIG. B A WHEEL SYSTEM WHOSE MOMENTUM CAN BE CONTROLLED WITH 
(IN THE LIMIT) ZERO ENERGY LOSS 
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APPENDIX B 


THE FEASIBILITY OF CONTROLLING MOMENTUM WITHOUT ENERGY CONSUMPTION 


The purpose of this appendix is simply to show that, ^ principle , 
the net momentum of a flywheel system can he varied as desired over an 
indefinite period of time with loss in energy which, in the limit, may 
he zero. To this end an all-mechanical, energy-conservative system is 
offered, for demonstration only, in Fig. B. 

Two coaocial reaction wheels are arranged with an idler-wheel drive 
between them. The idler is positioned hy a carriage, as shown in the top 
view. 


If the carriage position, x, is zero, then the two wheels run 
at equal speed, in opposite directions, and the net momenttua of the 
pair is zero. Let this initial wheel speed he ‘Xo). When x has a 
non-zero value the ratio of the reaction wheel speeds is given hy the 
geometry, and the net momentum of the system can have any value between 
(in the limit) plus and minus ‘J~2 JW(o), depending on x. 

For example, for the position shown, x = ^/P, the lower wheel must 
he running three times as fast as the upper wheel, and the net momentum 
of the system is 2JP-q_, directed downward. Moreover, if the system was 
moved from x = 0 to x = i/p without energy loss, then 
shown to be .45 ^(o) (and % = 1.34 P(o)), so that the net momentxmi 
of the system is .89 J<*^(o). 

When a change in x is made the necessary accompanying transfer 
of momentum, from wheels to vehicle, takes place, of course, via the 
torque impulse the idlers exert on the carriage (and thence on the 
vehicle) about the reaction wheel axis. (Unfortunately, the system 
shown also exerts a small torque impulse about an axis parallel to 
the idler axes. This must he eliminated, hy design, or compensated by 
the control computer, for three-axis operation. ) 

To show that operation with negligible energy loss is reasonable 
we consider, first, that there is viscous friction between the idler and 
the reaction wheels. Next we let the friction become high and note the 
limiting values of momentum change and energy change. Finally we consider 
the energy necessary to position the idlers. 

If the idlers in Fig. B are assumed to be massless, and if the 
friction between Idlers and reaction wheels is considered viscous, then 
the response of the system to a sudden change* in x, from xp to X 2 , 
turns out to be given by: 


* 

The sudden-change response is pursued because then x, which 
appears as a multiplier in the differential equations, is constant. 


- 44 


s + 


IT ~lfi) T 1")“ 

s[s + (l + Xg/ i^)o 


in which a is the reciprocal of the time constant of one reaction wheel 
skidding against a clamped idler located at radius i. 

If a large change in x is made suddenly from Xj_ to X2, then 
considerable energy is lost to friction. But if the change is made 
gradually, over a time which is long compared to the time constant of 
the idler system , l/o(l + x|/i^), then the energy lost goes to zero. 

That is, by making the momentum change gradually it can be made reversibly, 
in the limit. 

But since a is proportional to viscosity the speed with which 
momentum can be changed, with negligible energy loss, can be made as 
fast as desired by increasing viscosity as needed. 

We must consider also energy required to move the idlers. Here a 
possible refinement for reducing skidding is to make the idler wheels 
steerable, like the wheels of a car. Then a tiny force can steer the 
wheels rapidly to any desired position with energy loss which, again, 
can be made very small. 

The system presented in this appendix is not operationally optimum 
for several reasons. (For example, only l/2v2 times the total momentum 
capacity of the pair of wheels is utilized.) But the point is made, it 
is believed, that it is quite feasible to operate reaction wheel systems 
satisfactorily for long periods with the torquer system turned off. 
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